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Abstract 

An off-lattice Monte Carlo algorithm for solutions of equilibrium polymers 
(EP) is proposed. At low and moderate densities this is shown to reproduce 
faithfully the (static) properties found recently for flexible linear EP using a 
lattice model [Qj. The molecular weight distribution (MWD) is well described 
in the dilute limit by a Schultz-Zimm distribution and becomes purely expo- 
nential in the semi-dilute limit. 

Additionally, very concentrated molten systems are studied. The MWD 
remains a pure exponential in contrast to recent claims The mean chain 
mass is found to increase faster with density than in the semi-dilute regime due 
to additional entropic interactions generated by the dense packing of spheres. 

We also consider systems in which the formation of rings is allowed so 
that both the linear chains and the rings compete for the monomers. In 
agreement with earlier predictions the MWD of the rings reveals a strong sin- 
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gularity whereas the MWD of the coexisting hnear chains remains essentially 
unaffected. 

February 1, 2008 



PACS numbers: 82.35.+t, 61.25H, 64.60C 

I. INTRODUCTION 

The molecular mass distribution (MWD) of systems of linear unbranched "equilib- 
rium polymers" (EP) is essentially exponential In EP systems polymerization takes 
place under condition of chemical equilibrium between polymer chains and their respective 
monomers. A classical example we will focus on is provided by systems of surfactants form- 
ing polydisperse solutions of long worm-like aggregates, called "giant micelles" (GM), which 
combine with each other, or break into smaller fractions 0. 

Despite polydispersity, EP resemble in many aspects conventional quenched polymers 
where the polymerization reaction has been deliberately terminated. Recently, the basic 
scaling predictions for EP based on classical polymer physics have been tested by two 
of us (AM, JPW) III by means of a lattice Monte Carlo simulation based on the widely used 
Bond-Fluctuation Model (BFM) |^. This demonstrated excellent agreement with theory 
over a very broad range of density and temperature variation. Specifically, it was shown 
that the MWD takes the form, 

I exp{-x)dx (N* > (iV),0* < 0) 

c[x)ax oc < (Ij 

I x'^-^ exp(-7x)rfx (A^* < {N),(f)* > 0). 
The scaling parameter x = N/ (N) is the ratio of the chain mass A^ and the mean mass (A^), 
A^* and (p* mark the mean mass and the density at the crossover from dilute to semi-dilute 
regimes at given scission energy E, and 7 ~ 1.158 is the susceptibility exponent of the n —>■ 
vector model in 3D 0. The mean chain length (A^) was confirmed to vary with density 
and the (dimensionless) scission energy E as 
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{N) ^ N* (0/0*)" oc 0° exp{5E) (2) 

with exponents aj = Si = 1/(1+7) ~ 0.46 in the dilute and «// = (l + (7 — l)/(3z// — 1))/2 ^ 
0.6, 5// = 1/2 in the semi-dilute regime. The exponent z// ^ 0.588 is the swollen chain (self- 
avoiding walk) exponent in 3D. 

Recently, these results have been questioned in an interesting computational study where 
a reptation algorithm was used on a cubic lattice 0. It was suggested that the MWD 
becomes singular c{N) oc A^""^ with r 0.56 at very high volume fractions of order unity. 
This finding, if corroborated, might be of some relevance in view of the alleged Levy-flight 
dynamic behavior in system of GM 0. This claim prompted the present off-lattice Monte 
Carlo (OLMC) approach to be applied to the high density limit which was not accessible 
within our previous approach. We stress that effects only occurring at volume fractions of 
order one are unlikely to be relevant for real GM, but would be interesting on more general 
theoretical grounds in that the MWD in fact probes the free energy of the dense packing of 
beads. 

As we are going to show, new physics (i.e. additional terms in the system free energy) 
intervenes indeed due to packing effects of the spherical beads. This increases effectively 
the growth exponent a for (A^), but does not effect the scaling form of the MWD. We do 
not observe any trace of singularity — in perfect agreement with the analytical predictions 



The investigation of the high density limit is only one motivation for the off-lattice 
algorithm proposed. In addition to this an effective OLMC for EP algorithm is highly 
warranted to overcome the usual shortcomings of lattice models and to serve in examining 
the role of polymers semi-flexibility. It is also a better tool in dynamic studies of a broader 
class of soft condensed matter systems where bifunctionality of the chemical bonds might 
be extended to polyfunctional bonds, as this is the case in gels and membranes. Note that 
the OLMC was already applied successfully to rheological properties of EP reported in 
and to systems of EP brushes |T^ . 
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The scope of this paper is three-fold: we want to present the OLMC scheme and to test 
it by comparing the new results with previous findings (in the dilute and semi-dilute regime) 
obtained by means of our lattice Monte Carlo approach described in |Jl[|. Secondly, we wish 
to address the physics in the melt density regime. In addition, the OLMC algorithm is tested 
in systems in which the formation of closed loops is allowed. We demonstrate for the first 
time in a computer experiment that there is a singularity of the ring MWD co{N) oc A^~^ 
where r = 5/2 in 3D |ll],|l2|. 



This paper is organized as follows. After presenting the algorithm in Sec.|l| we focus in 



Sec.|T| on the different density and chain length regimes reflected, e.g., by the distribution 
of the radius of gyration {R'^{N)) versus mass N and other conformational properties. We 
discuss subsequently the MWD in Sec.|^. There we compare our OLMC results for mixed 
systems (in which rings are also present) with data obtained with the BFM and with a 
Grand Canonical lattice algorithm based on the mapping of the EP problem on a Potts 
model. The scaling of the mean chain mass is considered in Sec.[V|. The theoretical concepts 
(which have already been extensively considered elsewhere [§|,^) are briefly discussed en 
passant. We show that agreement of the OLMC with our previous work is excellent for 
small and moderate densities. For volume fractions larger 0.1 we evidence a third molten 
regime. In the final Section^ we summarize our findings. 

II. THE COMPUTATIONAL ALGORITHM 

It is clear that in a system of EP where scission and recombination of bonds constantly 
take place the particular scheme of bookkeeping should be no trivial matter ||13[. Since 



chains are only transient objects the data structure of the chains can only be based on the 
individual monomers, or, even better, on the saturated or unsaturated bonds This idea 
is depicted in Fig.|I](a). Each bond is considered as a pointer, originating at a given monomer 
and pointing to the respective other bond with which the couple forms a nearest neighbor, 
or to NIL (nowhere), if the bond is free (unsaturated). Chains consists of symmetrically 
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connected lists of bonds: jbond = pointer (ibond) ^ ihond = pointer {j bond). Recombina- 
tion of the two initially unsaturated bonds ibond = 2 and jbond = 5 connects the respective 
monomers imon = 2 and jmon = 5 in Fig.|I](a). Note that only two pointers have to be 
changed and that the remaining chains behind both monomers are not involved. Breaking a 
saturated bond ibond requires resetting the pointers of the two connected bonds ibond and 
jbond = pointer (ibond) to NIL. 

As mentioned in the Introduction this data structure has been incorporated ^ within 
the widely used BFM algorithm For the off-lattice version presented here we have now 
harnessed a very efficient bead-spring algorithm for polymer chains (for technical details 
see [|^) and cast it onto the data structure described above. This off-lattice Monte Carlo 
(OLMC) scheme is characterized by the bonded and the non-bonded interactions shown in 
Fig.|(b). 

Each bond is described by a shifted FENE potential where a bond of length r has a 
maximum at rmnT = 1 



UpENEir) = -K{rmax " ^o)^ lu 



r - ro 



E. (3) 



E corresponds to a constant scission energy. Note that UpENEij = ro) = —E and that 
Ufene near its minimum at ro is harmonic, with K being the spring constant, and the 
potential diverges logarithmically to infinity both when r — > r^ax and r — > r^in = 2ro — r^aa;- 
Thus the FENE potential does not need to be truncated at r^m and r^ax and discontinuities 



in the derivative of the potential are avoided. Following ref. we choose the parameters 
rmax — ro = ro — rmin = 0.3 and K/T = 40, T being the absolute temperature. The units 
are such that the Boltzmann's constant ks = 1. 

The non-bonded interaction between effective monomers is described by a Morse-type 
potential, r being the distance between the beads 

UM{r) = exp[-2a(r - rmin)] - 2 exp[-a(r - rmin)] (4) 

with parameters a = 24 and r^m = 0.8. The latter parameter sets roughly the sphere 



diameter. The 0-temperature of the coil-globule transition for our model is ~ 0.62 so 
that at T = 1 we work under good solvent conditions p^ . 

The model can be simulated fairly efficiently with a dynamic MC algorithm, as de- 
scribed previously |T^. The trial update involves choosing a monomeric unit at random and 
attempting to displace it randomly by displacements Ax, Ay, IS.z chosen uniformly from the 
interval —0.5 < Ax, Ay, A^; < 0.5. Moves are then accepted according to the Metropolis 
criterion and one Monte Carlo step (MCS) involves as many attempted moves as there are 
monomers in the system. 

In equilibrium polymers the bonds between neighbors along the backbone of a chain 
are constantly subject to scission and recombination events. In the present model only 
bonds, stretched a distance r beyond some threshold value, r^, are attempted to break so 
that eventually an energy UpENEi^) > could be released if the bond breaks. Since each 
monomer may have at most two bonds in the same time, all particles with unsaturated 
bonds (two for single monomers and one for chain ends) may form new bonds, once they 
approach each other within the same interval of distances Vf, < r < 1 where scissions take 
place. (Note that recombination for r < Vf, would violate detailed balance.) 

We do not allow in the presented study for branching. However, more than two bonds per 
monomer are possible in principle and this feature may readily included in the algorithmic 
framework. The generalization on netted structures like membranes or sponges is evident. 
(Note that this is less straightforward in the BFM scheme due to its lattice character which 
generates ergodicity problems.) Obviously, another big advantage of the off-lattice scheme 
compared to its lattice precursor is its applicability to rheological problems . 



In most parts of this paper we focus on systems where no formation of ring polymers is 
allowed. This condition has to be observed whenever an act of polymerization takes place. 
Because there is no direct chain information in the data structure this has to be done by 
working up the list of pointers (which adds only four lines to the source code). In physical 
time units the simulation becomes faster for higher E: the number of recombinations per 
unit time goes like exp{—E), but the chain mass only grows as exp(£'/2). Obviously, the 



algorithm becomes even faster for the mixed systems discussed in Sec. [IV B| where the ring 
closure constraint has been dropped. 



III. CONFORMATIONAL PROPERTIES 

The presented algorithm allows us to simulate a large number of particles at very modest 
expenses of operational memory. Most of the results in the present study involve 65536 
particles for number densities between cj) = 0.125 up to = 1.5. Note that our highest 
densities correspond to very concentrated solutions. This can be better seen from the volume 
fractions which vary between v(f) ~ 0.03 and 0.33. The latter value has to be compared with 
the (only slightly larger) hard-sphere freezing volume fraction ( "Alder transition" ) of about 
one half — it correponds, in fact, to a relatively dense hard-sphere liquid. Note also 



that dense globules of maximal 0^2 are formed by quenched monodisperse polymers in 



bad solvent 15 



The effective bead volume v ~ ttI^/Q ^ 0.22 used above was estimated with the measured 
mean bond length / ^ 0.75 from Tab.|I[ (Note that / is weakly decreasing with density 



IJI and that it is slightly smaller than the "bead diameter" given by the shifted Morse 
parameter rmm-) A similar, marginally larger volume of f ~ 0.25 may be obtained from the 
virial expansion for quenched polymers ||15 |. 



The scission energy E was varied in a similar range as in [||] from E = A up to E = 12 
to yield sufficiently strong chain mass variation. Note that the (N) remains always two 
orders of magnitudes smaller than the total particle number within the box. Hence, from 
our previous study |l| one expects finite box-size effects to be small. This was indeed born 
out by finite-size test performed by varying the box sizes from 16'^ over 32^ to 64^. Only the 
results from the largest boxes simulated for a configuration {E, (p) are reported here. 

Periodically the whole system is examined, the MWD c{N) (discussed in the subsequent 
section) and the distributions of the squared end-to-end distance {Rl{N)) and gyration 
radius, {Rg{N)) (averaged over all chains of a particular mass A^) are counted and stored. 
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(Periodic boundary conditions are implemented and interactions between monomers follow 
the minimum image convention. The computation of the conformational properties of the 
chains then imply a restoration of absolute monomer coordinates from the periodic ones for 
each repeating unit of the chain.) Moments obtained from these distributions are presented 
in Tab. I and Tab.y. 

The different simulational regimes are most easily demonstrated by the conformational 
changes with density and chain mass. We demonstrate that the conformational properties 
for EP follow the same universal functions as conventional quenched polymers. In Fig.|^ 
(i?2(Ar))i/2 and (i?^ (AT)) 1/2 are plotted versus A^. One configuration in the dilute regime 
and one in the concentrated limit (both at same E = 7) have been presented. In the first 
case one clearly sees the swollen coil exponent. From the upper dashed line the persistence 
length of swollen EP is estimated as 6 = Rg/N^'^^^ ^ 0.92. In the latter case only chains 
smaller than the excluded volume blob of size ^ size are swollen. Larger chains (A^ ^ g) show 
Gaussian behavior with R^./ ~ -R^ oc N^^'^. This slope was used to estimate the variation 
of the effective persistence length b = R^/N^'^ with cj). As can be seen from Tab.y b decreases 
with increasing density, hence with decreasing blob size ^ which we have estimated directly 
from the intercept of the two slopes at the radius of gyration. Note that ^ is relatively small 
for the densities computed. 

The averages over all chains of the mean-square end-to-end distance R^. and the radius 
of gyration Rg are plotted versus (A^) in Fig.|^. Swollen and Gaussian scaling behavior are 
again obtained for low and strong chain overlap respectively. Note that, e.g., the chains for 
= 1 are Gaussian for (A^) > 10. Thus the position of a particular configuration {E, 0) 
with respect to the crossover from dilute to semi-dilute regime can be determined. This is 
consistent with the value 0/0* given in Table | where 0* = (N) /{AitR^/3) as usual. 
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IV. MOLECULAR WEIGHT DISTRIBUTION 



We focus first on systems without rings, as in the rest of the paper, and consider then 
the effects due to the ring closure constraint by allowing linear chains and rings to compete. 

A. Systems without rings 

Fig. ^ displays a typical MWD c{N) obtained with the OLMC algorithm at high density 
(0 = 1.5). This is compared with BFM data. The c{N) is normalized such that = 
J2n ^c{N). Note that the free monomers are counted as chains of length = 1. Both 
curves display to high accuracy nice exponentials. This is a generic result for strongly 
overlapping chains — even at extremely high volume fractions. Notably, at variance with a 
recent finding no sign of singularity is observed. 

In order to corroborate this result we try to scale the MWD obtained for different densities 
and scission energies. In Fig.|^ we have plotted the (properly normalized) MWD versus the 
natural scaling variable x = N/ (N) . In the high density limit (main figure) data from 
two densities at various scission energies E (as indicated in the figure) collapses on a single 
'master distribution' c(x) = exp(— a;). In the dilute limit, i.e. for the non-overlapping chains 
shown in the inset, we find again a data collapse, but with slightly different slope in linear-log 
coordinates. (Obviously, for large chain length the statistics deteriorates.) 

How can these results be rationalized? Within a standard Flory-Huggins mean-field 
approach one may write the total free energy density as 

oo 

^[cm = J2 <N) (log(c(iV)) + /iiV + UUN, 0, E)) . (5) 

N=l 

The first term on the right is the usual translational entropy. The second term entails a 
Lagrange multiplier which fixes the total monomer density (p. The most crucial last term 
encodes the free energy of a reference chain of length N in the field created by the surrounding 
chains and free monomers, f chain will in general depend on the chain length A^, the density 
and the bonded and non-bonded interaction parameters of the model studied . We have 
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not computed here irrelevant additive terms (such as virial terms) which are not conjugated 
to c{N). By functional derivation with respect to c{N) one readily obtains the equilibrium 
MWD: 



c(iV) = exp(-^ - Ud{N, 0, E) - fiN). 



(6) 



We have defined here fend = f chain — E — 1 and absorbed all contributions to fend linear in 
within the Lagrange multiplier. Hence, within the mean- field approximation eq.(^ the 

MWD discussed above probes directly the free energy of an EP chain (besides the terms 

linear in mentioned above which are fixed by the imposed density) which essentially 

"renormalizes" the scission energy E. 

We infer from Fig.^ and Fig.|^ that for our EP chains fend is to very high accuracy mass 

independent. Two classical examples where this is indeed rigorously true are non- interacting 

rigid rods and Gaussian chains 0. (There fend is even independent of 0.) If this is the case, 

the MWD becomes then a pure exponential 



where the scaling variable is x = N/{N), i.e. the inverse Lagrange multiplier equals the 
mean-chain length, {N)fi = 1. 

This is in agreement with our scaling in the high density limit which implies that fend 
is within numerical accuracy independent of A^. We have explicitly checked that here the 
MWD scales like c{N) oc exp{—E — fend{.<P) ~ ^1 {^)) where fend only depends on 0. The 
slightly different slope found in the dilute regime suggests, however, a weak mass effect. It is 
relatively simple to understand these results from the standard theory of polymers in good 
solvent solutions There the chemical potential of ends is given by fend = (7 — l)/'^/ log(0 
where ^ is size of the 'blob', i.e. the excluded volume correlation length for chains of length 
at density 0. 

In the semi-dilute limit we have ^ oc g^^ oc (^g denotes the number of 

monomers within a blob), hence 



c{N)dN/(p = c{x)dx = exp{—x)dx 



(7) 
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fend = (7 - l)/(3/// - 1) log(0/0o) (8) 

is independent of mass and scission energy. The constant reference density 0o was introduced 
here for dimensional reasons. Obviously, eq. (§) can strictly hold only in the asymptotic limit 
of large blob sizes (0-^0). Additionally, one can not exclude weak chain- length dependence 
of fend &t N ^ 1, but this becomes irrelevant in the scaling limit of (A^) ^ 1 |]18[ . 

In the opposite dilute limit the correlation length becomes density independent and is 
given by the size of the chain ^ = i? oc A^''. Hence, fend = (7 ~" 1) log(^) the MWD 
becomes the Schultz-Zimm form given in the Introduction (eq.|I]). Note that the 7-exponent 
in 3D is only slightly larger than its mean field value 7 = 1. Hence, the predicted power law 
depletion in the MWD for small A^ is very weak and requires a relatively large mean chain 
mass (A^) ^ 1 This is why within the range of the parameters accessible essentially only the 
slightly different exponential tail c{x) oc exp(— 7x) is visible. 

Obviously, some of our simulations are in the intermediate regime (0 ~ 0*) between 
these limiting cases (depicted by the two slopes in Fig.^. In order to characterize all the 
MWD with one parameter we fitted the tail of log(c(A^)) with —E — fi — fiN where both n 
and /i are fit parameters. In the dense limit naturally one finds yu(A^) = 7e// ~ 1. As soon 
as the overlap decreases 7e// —>■ 7. Note that the /i-variation with E for given is weak and 
we have only tabulated the value for the highest E, i.e. for the most strongly overlapping 
chains. The values are given in Tab.|l| and are discussed in the subsequent Section]^. 

B. Mixed systems with linear chains and rings 

First we briefly assess the importance of the ring closure constraint. We allow the 
formation of rings so that both linear chains and rings may coexist and compete for the 
monomers. In Fig.^ we present the MWD of rings Co(A^) and of linear chains Ci{N) at 
density = 1.5 and E = 7. The linear chains appear to be unaffected by the presence of 
rings and scale, as before, Ci{N) oc exp(— /lA^) where the slope /x is given by the inverse of 
the mean mass of the linear chains A^i. This can be easily understood by generalizing the 
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Flory-Huggins expression eq.(^ as a sum over the two different species s = 0, 1. The two 
different MWD decouple and the functional derivation with respect to both yields simply 
Cs{N) = exp(— /4aj„(A^, 4>, E) —Nfi) for both species in analogy with eq.(^. The free energy 
of the linear chain is f chain = E + 1 + fend as before. In the strong overlap regime fend 
becomes again independent of mass ( hence, the pure exponential seen for ci{N)) and, 
as a matter of fact, equals the density dependent value /i(0) obtained for the purely linear 
systems which is tabulated in Tab.|T|. Therefore the mean mass of linear chains determines 
once again the chemical potential fi. 

However, our systems of flexible rings appear to be ring dominated, i.e. much more 
monomers are contained in closed loops than in linear chains, as expected. In contrast 
in systems of semi-flexible polymer chains the formation of rings is going to be strongly 
suppressed [0]. The MWD of rings is not exponential, but rather shows a clearly pronounced 



power law behavior co{N) oc A^""^ where r = 2.5. (This does not exclude, however, a final 
exponential cut-off, which can hardly be detected because of the dominating power law.) 
The exponent r can be obtained by a simple argument due to Porte |TT]]. The ratio of both 
MWD gives the ratio of the two partition functions 

^ = ^ = exp(fi -f ) (9) 

^ (7\r) Z (iV) t'KJ chain J chain/' \ J 

The ratio Zi/Zq is equal to the probability of opening a loop, which must be proportional 
to: (i) the Boltzmann weight, exp (—£'), due to the constant scission energy, (ii) the number 
of places where the ring can break, A^, and (iii) the volume oc N^'^ that two neighboring 
segments can explore after being disconnected {D denotes the dimension and u the relevant 
Flory exponent at density (p). Hence, Cq oc N^'^ exp{~N/Ni) with r = Du + 1 = 5/2 at 
strong chain overlap in D = 3, and t = 2 in D = 2 where u = 1/2. Both exponents are well 
born out by the measured Co-slopes depicted in the inset of Fig. ||. We compare here data we 
obtained with three different algorithms: OLMC, BFM and Grand Canonical Potts model 
(in 2D and 3D). The exponents in both dimensions are fully consistent with the predicted 
ones. 
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We return now again to systems containing linear chains only. 

V. THE SCALING OF THE AVERAGE CHAIN MASS 

So far we have reported on the general form and scaling of the MWD. Now we want to 
go further and to investigate the scaling of the average chain mass and the second moment, 
i.e. the polydispersity index / = {N"^) / (N)"^ , with regard to density and scission energy 
E (see Tab.l). 

In the high density limit, i.e. within the validity of eq.(0), the inverse Lagrange multiplier 
equals the mean-chain length as mentioned above and the polydispersity 1 = 2. From the 
normalization constraint (p = J2n Nc{N) one infers that the mean mass is generally given 
by 

(AT) = = y/0exp(E + /,„,(0)). (10) 

For semi-dilute polymer chains one obtains from eq.(0) and eq.(|D (A^) oc 
^0i+(7-i)/(3^'/-i) exp(ii^) in agreement with the exponents an ~ 0.6 and 5[[ = 1/2 quoted 
in the Introduction. Similarly, one finds in the dilute limit «/ = 5/ = 1/(1 + 7) ~ 0.46 

In Fig.0 we have plotted (A^) versus E to check for the expected exp{SE) behavior. 
Despite the small difference in both values one distinguishes for dilute systems a slope 
with 61 = 0.46 and for the dense system (0 = 1) 611 = 0.5. One can verify on Fig.|^ 
that at concentration (p = 0.125 a crossover into the semi-dilute regime for E ^ 9 occurs 
whereby the slope of the exponential function changes from 61 ~ 0.46 to 6n = 1/2, i.e. the 
isolated polymer coils get large enough so they start touching each other. At = 0.25 the 
concentration is already sufficiently high so this happens at comparatively lower energies 
E > 5. The 6jj = 1/2-slopes confirm that at high densities fend becomes chain length- 
independent and that eq.(|iy) holds. This scaling prediction is verified explicitly in Fig.|^ 
where we have plotted (A^) versus the scaling variable indicated by eq. (p^ . Here we have 
used the values /i(0, E) determined directly from the MWD. A plot using the prediction for 
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semi-dilute polymers eq.(P) is, however, not successful for > 0.5. Hence, the measured /i 
encapsulates physics other than the one expected in the limit of large blobs. To corroborate 
this further we plot in Fig.^ u = (N)'^ / {(p exp{E)) versus (p. From eq.(|TO|) we know that 
u = exp(/e„d) for strong overlap. Not surprisingly the data collapse in that limit, but not 
the dilute systems. (This will be improved in Fig.|TD|, see below.) Also indicated are the /i 
(stars) directly estimated from the MWD in the strong chain overlap limit. Notably, for 
> 0.5 the mean mass appears to increase much faster than the growth exponent an = 0.6. 
This breakdown of eq.(||) — but not of eq.([lOl) — is not unexpected if one bears in mind 



that the scaling arguments can only be valid in the limit of — ^> and g — > oo, i.e. 
when the blob is larger enough. This result suggests to rewrite eq.(H) as a systematic series 
expansion 

fendicp) = 5_i log(0/0o) + 5i0 + ... (11) 

where in order to match semi-dilute and melt regime we choose B_i = (7 — 1)/ {Suj — 1) ^2.1 
and 00 ~ 0.018. The fit using eq.(|TT]) with Bi ^ 0.8, depicted in Fig.|^, is only qualitatively 
satisfactory and further study is warranted. The above expansion is motivated by a recent 
second virial theory on wormlike micelles and rigid rods interacting via hard-core excluded 
volume P]. The applicability of this theory to our systems of flexible EP with persistence 
length b of order of the sphere diameter (see Tab.|ID appears to be unclear. Nevertheless, it 
is interesting to note that the second virial coefficient Bi = 8/3 v ~ 0.6 predicted in is 
relatively close to the constant measured. Alternatively, the dependence of (A^) on might 
be represented by an effective growth exponent Oe// ~ 1 (dotted line). These values agree 
also favorably with an earlier lattice simulation although this is presumably accidental 
due to the different microscopic physics which must intervene. 

A proper crossover scaling between dilute and semi-dilute regimes (but not with regard 
to the molten regime) is achieved by plotting {N)/N* versus 0/0* in Fig.p!0|. The crossover 
length A^* oc exp{—ipE) with ip = (1/2 — ai)/{ajj — aj) ~ 0.26 and the crossover density 
0* oc exp{K,E) with k, = ai{l — ip) 0.34 are readily found by matching the asymptotic 
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behaviors of the dilute and semi-dilute regime. We have included data from both OLMC 
and the BFM. For clarity, the BFM data have been shifted downwards by a factor 10. As 
one anticipates from Fig.|^, only OLMC data for cj) < 0.5, i.e. volume fractions smaller 
0.1 collapse properly on the predicted asymptotes — confirming hence the polymer physics 
expressed in eq.(^) — while at higher densities simple liquids physics becomes relevant. 

VI. SUMMARY 

We have proposed here a new off-lattice Monte Carlo algorithm (OLMC) for systems 
of equilibrium polymers (EP). It is shown that this model faithfully reproduces the results 
of ref. |l|: the MWD of linear EP is essentially exponential even in the limit of very high 
density in contrast to recent claim 0]. Note that our findings do not support recent claims 
of Levy-flight dynamic behavior in EP and CM which would require a singularity in the 
MWD of linear chains 0] . 

If ring formation is allowed, however, the MWD for rings alone is strongly singular, 
cq{N) oc N^'^ with T = Du + 1. This result has been confirmed both in two- and in three 
dimensions for different lattice and off-lattice models. We have also shown that the MWD 
of linear chains is not affected by the presence of rings in the system which demonstrates 
that the general Flory-Huggins approach is very accurate even in this case. 

The mean chain length of linear EP is found to vary with density and scission energy 
E like (A^) oc </)"exp(5i?) as observed in our earlier Bond Fluctuation Model investigation 
lU with exponents a/ = 5/ ~ 0.46 in the dilute regime and an = 0.6, 6ii = 0.5 in the 
semi-dilute regime. This holds only for volume fractions smaller 0.1. 

At higher densities when the blobs become too small the scaling approach breaks down. 
Conformational properties are then largely determined by packing effects as in simple liquids. 
The mean chain length then grows much faster with 0, qualitatively following a nonalgebraic 
dependence similar to a recent decription put forward for rod-like and semi-flexible micelles 
Our result is also compatible to an effective growth exponent cte// ~ 1, in agreement 
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with recent simulations and even with experimental observations although in the 
latter case this agreement might be accidental in view of the rather low volume fraction 
probed experimentally. 

As a further development the OLMC algorithm the model will be used for studies of the 
dynamic behavior of linear and ring EP for which investigations are currently underway. 
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TABLES 





E 


{N) 


I 


Re 


Rg 




0.125 


8 


25 


1.86 


6.4 


2.5 


0.3 


0.125 


9 


42 


1.88 


8.7 


3.4 


0.5 


0.125 


10 


66 


1.89 


12.2 


4.5 


0.7 


0.125 


11 


108 


1.93 


14.8 


5.9 


1.0 


0.125 


12 


181 


1.95 


19.9 


8.0 


1.5 


0.25 


7 


24 


1.88 


6.0 


2.4 


0.6 


0.25 


8 


38 


1.90 


7.9 


3.2 


0.9 


0.25 


9 


63 


1.93 


10.4 


4.2 


1.2 


0.25 


10 


102 


1.93 


13.6 


5.5 


1.7 


0.25 


11 


166 


1.92 


17.5 


7.1 


2.3 


0.25 


12 


270 


1.99 


23.1 


9.3 


3.1 


0.5 


7 


38.4 


1.92 


7.9 


3.2 


1.8 


0.5 


8 


62 


1.94 


10.3 


4.2 


2.5 


1 


5 


28 


1.92 


5.8 


2.3 


1.9 


1 


6 


46 


1.94 


7.5 


3.0 


2.5 


1 


7 


74 


1.94 


9.6 


3.9 


3.4 


1 


8 


123 


1.96 


12.7 


5.2 


4.8 


1.38 


7 


105 


1.97 


10.9 


4.5 


5.0 


1.5 


6 


73 


1.95 


8.6 


3.6 


4.0 


1.5 


7 


120 


1.97 


10.6 


4.6 


5.1 



TABLE I. Summary of measured quantities for configurations with (A) > 20. Quantities 
tabulated: the mean chain mass (A), the polydispersity index / = (A^)/(A)^, the mean end-to-end 
distance Re the mean gyration radius Rg, and the chain overlap <^/0*. Note that / increases 
systematically with chain overlap. 
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b 


s 


J i 


0.125 


0.7574 


1.41 


4.7 


0.4 


0.25 


0.7565 


1.33 


4.0 


0.65 


0.5 


0.7530 


1.2 


3.8 


0.9 


1.0 


0.7408 


1.12 


3.3 


1.66 


1.38 


0.7234 


1.05 


2.8 


1.97 


1.5 


0.7158 


1.02 


2.5 


2.21 



TABLE II. Density dependence of mean bond length / = (/^)^^^, the persistence length b for 
Gaussian chains of blobs, the blob size £ estimated from the slopes in Fig.^, and the free energy 



factor /i((/>) obtained from the exponential tail of the MWD as explained in Sec. IV. 
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FIGURES 




FIG. 1. Sketch of algorithm: (a) Bonds between spherical beads break and recombine constantly 
with rates depending on the scission energy E which is assumed independent of mass N and 
number density (f). Each monomer has two (saturated or unsaturated) bonds. Chains consists of 
symmetrically connected lists of bonds. The data structure is based on the bonds rather than the 
polymer chains, (b) Plots of bonded (FENE) and non-bonded (Morse) interactions used in the 
present model. The shaded area denotes distances where scission - recombination events may take 
place. 
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r v=0.588 

10° 10' 10' 10^ 

N 

FIG. 2. Variation of the end-to-end distance Re (triangles) and the radius of gyration Rg 
(squares) with chain mass N at E = 7 for dilute (open symbols: E = 7,(p = 0.125) and concentrated 
(full symbols: E = 7,(j) = 1) systems. The dashed lines denote the dilute exponent = 0.588 
which is visible for the dilute systems over the full range of and for the semi-dilute system for 
small A^, i.e. within the blob. The slope i^u = 1/2 indicates the Gaussian statistics for strongly 
overlapping long chains. 
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FIG. 3. End-to-end distance Re and radius of gyration Rg vs. mean chain mass (A^). At small 
density ((/> = 0.125,0.25) and, hence, weak chain overlap the chains are swollen (dashed line). At 
high density and mass we obtain Gaussian statistics (solid lines). 
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FIG. 4. MWD c(A^) versus at high density for OLMC (upper curve; E = 7,(1) = 1.5) 
and BFM (lower curve; E = 7,8(p = 0.5) systems. Note that the distribution is always a pure 
exponential and no sign of singularity was found for whatever density or energy. These curve are 
used to fit /i. For OLMC we have (N) ^ 119, fi w 2.21, for BFM (N) w 98, fi w 2.85. 
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x=N/<N> 

FIG. 5. Data collapse of (properly normalized) MWD c{x) versus x = N/{N) for OLMC. In 
the main figure the high density scaling prediction c(x) = exp(— x) is verified (full line). The open 
symbols denote density (f> = 1 and the full symbols 4> = 1.5. Inset: The MWD in the dilute limit 
(OLMC at (/) = 0.125) compares well with the prediction c(x) oc exp(— 7x) (dashed line). The same 
symbols are used as in the main figure. 
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N 

FIG. 6. MWD for mixed systems containing rings. Main figm'e: cq-MWD of rings (spheres), 
ci-MWD of linear chains (squares). Inset: Ring MWDs for different models and dimensions. The 
indicated slopes confirm r = 5/2 in 3D and r = 2 in 2D. 
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FIG. 7. Variation of mean chain mass (A^) with dimensionless bond energy E for various 
number densities (/> as indicated in the figure. At low chain overlap the data are consistent with 
the exponent 5/ = 0.46 (dashed lines), at hight density with the mean field exponent bn = 0.5 
(full lines). 
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FIG. 8. Scaling attempt in the high density hmit for (A^). The coUapse confirms the scahng 
with respect to u in the high chain overlap limit and the f\ obtained independently from the MWD 
c{N). Note that the remarkable good scaling works even in the dilute limit. The dashed line marks 
for comparison the dilute exponent a/. 
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FIG. 9. A crude rescaling attempt of mean chain mass u = (A^)^/((/)exp(i?)) versus (f). Also 
indicated are the exp(/i) estimated from the MWD (stars). The mean chain mass rises faster for 
large densities (p ^ 0.5. This regime is compared with an effective exponent Oe// ~ 1 (dotted line) 
and the nonalgebraic density dependence eq.(0) with B^i = 0.21 and Bi ~ 0.8 (dashed-dotted 
line). 
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FIG. 10. Dilute-Semidilute crossover scaling {N)/N* vs. for OLMC and BFM (shifted 

downwards for clarity). For densities (f) > 0.5 the OLMC data do not scale and are systematically 
above the semi-dilute asymptote. This is due to additional entropic interactions between the dense 
beads which increases the mean chain mass. 
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